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ABSTRACT 

We review three main results of our recent study: 

• We show that a proper treatment of the tidal interaction prior to the onset 
of the common envelope (CE) leads to an enhance mass loss. This might 
increase the survivability of planets and brown dwarfs that enter a CE phase. 

• From the distribution of planets around main sequence stars, we conclude 
that around many sdB/sdO stars more than one planet might be present. 
One of these might have a close orbit and the others at about orbital periods 
of years or more. 

• We show that the intense ionizing flux of the extreme horizontal branch star 
might evaporate large quantities of a very close surviving substellar object. 
Balmer emission lines from the evaporated gas can be detected via their 
Doppler shifts. 



1. Introduction 

EHB (Extreme Horizontal Branch) stars are hot, small, helium burning stars. In order 
to become EHB the RGB progenitor must lose most of its envelope. For the purpose of these 
proceedings there will be no differentiation between sdB, sdO (which are the spectroscopical 
classes) and EHB stars (the photometric definition). It is known that planets can exist 
around sdB stars. In many cases it is possible that planets are responsible for the formation 
of the EHB star, e.g. HD 149382 b (Geier et al. 2009). In these proceeding we will review 
three main topics: In section 1, we will discuss the interaction of brown dwarfs or low main 
sequence stars and the RGB progenitor. Our goal is to study in more detail the evolution 
of binary systems in a stage prior to the onset of the CE phase, and in particular systems 
that have reached synchronization; the synchronization is between the orbital period and the 
primary rotation period. In section 2 we will discuss the bimodal distribution of planets, and 
the implications to EHB stars. In section 3 we will discuss planet evaporation, and detection 
of Ha, and H/3 emission. In section 4 we will summarize our concussions. 
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2. Interaction between brown dwarfs or low main sequence stars and the 

RGB progenitor 

We start our calculation when the primary stellar radius has increased enough for tidal 
interaction to become significant. For the binary systems we study, where the primary is an 
RGB star and the secondary is a low-mass main sequence (MS) star or a brown dwarf, tidal 
interaction becomes important when the giant swells to a radius of Rg ~ 0.2a, where a is the 
orbital separation and Rg is the giant radius (Soker 1998). The primary radius increases along 
the RGB as the core mass increases. The primary and the secondary initial masses range 
in the binary systems we study are O.SMq < Mi < 2.2Mq, and O.OISMq < Ms < 0.2Mq, 
respectively. 

When a binary system starts its evolution it is not synchronized, and therefore tidal 
interaction will lead to a fast spiraling-in process. The binary system can then either reach a 
synchronization or stays asynchronous. In systems that maintain synchronization two oppos- 
ing forces act: On one hand in order to maintain synchronization the secondary transfers, via 
tidal forces, orbital angular momentum to the envelope, and the orbit shrinks. On the other 
hand mass loss acts to increase orbital separation. If orbital separation increases faster than 
the RGB stellar radius, no Common Envelope (CE) will occur either due to total envelope 
loss or to a Helium fiash. However, if orbital separation decreases relative to the RGB stellar 
radius, and He fiash or total envelope loss do not occur too early the secondary enters the 
envelope either due to Darwin instability, or by the swelling RGB envelope. By that time 
the envelope mass is lower the the initial envelope mass. These processes are represented in 
figure [H 

As can be seen from the figure tidal interaction before the formation of the CE increases 
the likelihood of low mass companions to survive the CE phase. Furthermore, higher mass 
loss rate decreases the chance of a late CE formation, but if late CE occurs, it does so with 
lower envelope mass (for details see Bear & Soker 2010a). 

3. Bimodality of planets around MS stars 

We study the distribution of exoplanets around main sequence (MS) stars and apply our 
results to the binary model for the formation of extreme horizontal branch (EHB; sdO; sdB; 
hot subdwarfs) stars. By the binary model we refer not only to stellar companions to RGB 
stars (Han et al. 2002; Han et al. 2003), but to substellar objects as well (Soker 1998). The 
bimodal distribution presented in figure |2] is taken from the Planets Encyclopedia edited by 
Jean Schneider; http : // exoplanet .eu / ) . 
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Normal mass loss rate 




M2{Mq) 

Fig. 1. — Mg vs. M2 representing binary systems that reached or did not reach synchroniza- 
tion. Our normal mass-loss rate oi rjn = 3 x 10^^ was used (for details see Reimers 1975, & 
Bear & Soker 2010a. The initial core mass is Mc(0) = OAMq, and the initial (primordial) or- 
bital separation prior to tidal interaction is Uq = 5Rg{0) = 405Rq. Calculation is terminated 
as indicated, when one of the following occurs: The Darwin instability brings the system to 
a CE phase (marked Darwin Instability); The core mass reaches O.48M0 and assumed to go 
through a core Helium flash (Helium flash). The two other channels: (1) A total depletion 
of the RGB envelope. (2) The RGB stellar radius exceeds the orbital separation {Rg = a), 
do not occur for the case presented here. 
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Fig. 2. — Number of planets vs. log 
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In this work we follow Soker & Hershenhorn (2007). Soker & Hershenhorn (2007) 
examined the number of planets as a function of metallicity bins and the planet mass Mp, 
orbital separation a, and orbital eccentricity e, in several combinations. They found that 
planets orbiting high metallicity stars tend to part into two groups in a more distinct way 
than planets orbiting low metallicity stars. They also found that high metallicity stars tend 
on average to harbor closer planets. Soker & Hershenhorn (2007) had 207 planets in their 
analysis. We repeated their analysis using 331 planets (out of more than 400 that were 
discovered so far) and got similar results. In figure [2] only the multi-planet systems from the 
research are represented (for more details see Bear & Soker 2010b). 

The planets shown in this figure are all part of a multi-planet systems (two planets or 
more that orbit the same star). All central stars are in the mass range of 1 — 1.5Mq, this 
mass is the typical mass of the progenitors of EHB stars (for more details see Bear & Soker 
2010b). In general, there are two groups of planets, the inner marked by (I), and the outer 
planets that are marked by (O). In order to implement this distribution to EHB stars we 
will review the observations of planets around EHB stars: 

Geier et al. (2009) announced recently the discovery of a close substellar companion 
to the hot subdwarf (EHB) star HD 149382. The orbital period is very short, 2.391 days, 
implying that the substellar companion had evolved inside the bloated envelope of the pro- 
genitor RGB star (a CE phase). The mass of the companion is 8 — 23Mj , so either it is a 
planet or a low mass brown dwarf. This discovery supports the prediction of Soker (1998) 
that such planets can survive the common envelope (CE) phase, and more relevant to us, 
that planets can enhance the mass loss rate on the RGB and lead to the formation of EHB. 

However, Jacobs et al. (2010) analyzed He lines and found no evidence for the presence 
of this claimed planets. This debate will soon be resolved by further observations. Other 
planets that orbit EHB at larger separations have been detected (Silvotti et al. 2007; Lee 
et al. 2009 & Qian et al. 2009). Silvotti et al. (2007) announced the detection of a planet 
with a mass of 3.2Mj and an orbital separation of 1.7AU, around the hot subdwarf V391 
Pegasi. Serendipitous discoveries of two substellar companions around the eclipsing sdB 
binary HW Vir at distances of 3.6AU and 5.3AU (Lee et al. 2009) and one brown dwarf 
around the similar system HS 0705-1-6700 with a period of 2610(i and a separation of < 3.6AU 
(Qian et al. 2009). Recently Geier et al. (2010) discovered a brown dwarf companion to 
the hot subdwarf SDSS J083053. 53-1-0000843. 4. This system contains an sdB star with an 
approximated mass of 0.25 — O.47M0. A brown dwarf of 0.045 — O.O67M0 orbits this sdB 
with an orbital period of 0.096d. Due to its close orbit it is very likely that this system went 
through a Common Envelope phase. 

All of these five systems are present in the graph, with their evaluated orbital separation 
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around the progenitor of the EHB star (Bear & Soker 2010b). It is quite plausible that closer 
planets did interact with the RGB progenitor of the sdB star; they are not observed in these 
systems. We end by noting that all these substellar companions have been detected in the 
field. The main conclusions we can draw from figure [3] are as follows: Planets are expected 
in a double peak distribution. Outer planets can survive the evolution. In particular, inner 
planets in the same system that were engulfed, saved the outer planets by enhancing mass 
loss rate early on the RGB. The inner planet might survive the CE phase and be found 
around the EHB star, but only if massive enough M > lOMj. 



4. Planet evaporation and detection 

We study the evaporation of planets orbiting EHB stars. We adopt the simple model 
presented by Lecavelier des Etangs (2007) which represents the blow-off mechanism (Erkaev 
et al. 2007) and investigate the implications for a planet orbiting an HB star (this model is 
similar to the energy limited model purposed by Murray-Clay et al. (2009). We refer to the 
ionization model as well as a lower limit for the mass loss from the planet (for details see 
McCray & Lin 1994). 

When the central source is hot a large fraction of the radiation is energetic enough to 
ionize the evaporated gas. The evaporated gas recombines and emits at longer wavelength, a 
radiation that escapes from the planet's vicinity. Although recombination is not relevant to 
planet around solar-like stars, its role becomes more important for hot HB stars and central 
stars of planetary nebulae. We assume that: 

• Most of the evaporated gas flows toward the radiation source. 

• The central star keeps the gas almost fully ionized. 

• The ionizing photons of the parent star that are absorbed by the evaporated gas are 
removed from the radiation that heat the star. 

• All the radiation emitted by the recombining gas escape. 

• We assume that the gas flows with the escape velocity from the planet. 

The recombination rate is proportional to the density square, hence to mass loss square. 
We solve the mass loss rate of Lecavelier des Etangs (2007) taking into account the recom- 
bination. Substituting numerical values gives 
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Hp [AU] 

Fig. 4. — Mass evaporation rate rhp (left axis) versus the orbital separation ap. The right axis 
gives the total mass that would be evaporate during a period of 6 x lO^yr. The blue circles 
(lower line) represent the ionization model (for details see Bear & Soker 2010c). The black 
thick (upper) line represents the evaporation rate based on Lecavelier des Etangs (2007) for 
a black body energy distribution. The red thick line represents the evaporation rate based on 
Lecavelier des Etangs (2007) for a self consistent calculated spectrum of HD 149382 (Heber 
2010). The blue thin line represents the same model of Lecavelier des Etangs (2007) with 
recombination of the evaporated gas included (eq. [T]) for a self consistent calculated spectrum 
of HD 149382 (Heber 2010), instead of a black body that is not accurate at wavelengths below 
1200A. The evaporation rates are calculated for an EHB central star and a planet with the 
properties of the HD 149382 system: Tehb = 35500ir, Mehb = O.SM©, i?EHB = O.URq, 
Mp = 15Mj (Geier et al. 2009) and Rp = O.IRq. The orbital separation of this system 
is Op = 5 — 6.1i?0, but here it is an independent variable. The magenta line represents an 
orbital separation of Up = 5.5Rq. 
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Where An/S is the sohd angle to where the evaporation flow occurs with /3 = 0.5, v^sc is 
the escape velocity of the gas, Rp is the planet radius and ~ 20eV is the average energy of 
the ionizing photons. Figure H] represents the different mass loss considered and the mass loss 
that takes into account recombination. These are calculated with the appropriate spectrum 
as was calculated for HD 149382 (Heber 2010). For comparison we show the evaporation 
rate for a BB (Black Body) spectrum with the same effective temperature and luminosity 
(black upper line, for details see figure caption). 

The properties of the EHB central star and the planet are taken to be those of the HD 
149382 system (Geier et al. 2009; see figure caption). The orbital separation of this system 
is Op = 5 — G.IRq, but in the figure this is an independent variable. On the right axis of 
figure m we give the total mass that would be evaporate during a period of 6 x lO^yr, about 
the duration of the HB phase, with the same mass loss rate given on the left axis. We will 
concentrate on the orbital separation range of Up = 0.01 — 5AU. 

The calculation of the Ha luminosity from the evaporated gas is done in the following 
way, starting with the following assumptions: 

1. The evaporation is mainly into a solid angle 47r/3. 

2. Close to the planet, where most of the recombination occurs, the material flows at the 
escape speed from the planet. 

3. For typical values we find the medium to be optically thin to Ha. 

4. We assume that the evaporated gas is almost completely ionized. Any recombination 
that occurs is balanced by the incoming photons from the EHB star. 

5. Most of the recombination and the Ha source occur at a relatively high density of n ~ 
10^*^ — lO^^cm"^. At such densities collision between atoms change the amount of energy 
that is channelled to Ha- In our simple treatment we take the recombination coefficient 
neglecting the dependence on density. We note that Bhatt (1985) calculates the Ha 
emission from a destructed comet. He estimates the density to be ~ lO^'^cm"^ and 
neglects the dependence on density. Korista et al. (1997) found that the dependence 
in density on this high densities is negligible. 

The Ha energy released due to recombination is 




(2) 
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Solving the integral yields 

We find equivalent width of EWa ~ 0.05A for the Ha emission and EW^ ~ 0.006A for 
the Hf5 emission, both for the calculated spectrum of (Heber 2010). 

Although the EWs are not high, their periodic variation might ease the detection of the 
line. At an orbital separation of 5.5Rq the orbital velocity of the substcllar companion is 
130km s~^. Therefore, during the orbital period the center of the emission by the evaporated 
gas might move back and forth over a range of up to ~ 5.5A and ~ 4. OA, for the Ha and 
Hl3 emission lines, respectively. These EWs are an upper value since they are based on the 
black body distribution. We conclude that it might be possible to identify a planet via the 
Ho; emission of its ablated envelope. 

5. Summary and Conclusions 

We find that the pre CE evolution is crucial in understanding binary systems where 
the primary star is an evolved red giant branch (RGB) star, while the secondary star is a 
low-mass main sequence (MS) star or a brown dwarf. Moving to planets, from studying the 
distribution of planets we see that they are expected in a double peak distribution. Outer 
planets can survive the evolution of the progenitor of the EHB. In particular, inner planets 
that were engulfed by the RGB progenitor might "saved" the outer planets by enhancing 
mass loss rate early on the RGB. With the enhanced mass loss rate the RGB star will form 
an EHB star. We also saw in section 3 that planets close to an EHB star might be detected 
through Ha and H/3 emission. 
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